Chapter 1 

(N 



Supersymmetry in Random 
Matrix Theory 



S 

PL,- 

-I— > 

S: 

Thomas Guhr 1 

1 Fakultat fiir Physik, Universitat Duisburg-Essen, 
r — | Lotharstrasse 1, 47057 Duisburg, Germany 

Q\' 

. ■ Abstract 

in 

o, 

Supersymmetry is nowadays indispensable for many problems in Random Ma- 
trix Theory. It is presented here with an emphasis on conceptual and struc- 
tural issues. An introduction to supermathematics is given. The Hubbard- 
Stratonovich transformation as well as its generalization and superbosonization 
are explained. The supersymmetric non-linear a model, Brownian motion in 
superspace and the color-flavor transformation are discussed. 

1.1 Generating Functions 

We consider N x N matrices H in the three symmetry classes [Dys62] real 
symmetric, Hermitean or quaternion real, that is, self-dual Hermitean. The 
Dyson index f3 takes the values f3 = 1,2,4, respectively. For j3 = 4, the iV x 
iV matrix H has 2x2 quaternion entries and all its eigenvalues are doubly 
degenerate. For a given symmetry, an ensemble of random matrices is specified 
by choosing a probability density function P{H) of the matrix H. The ensemble 
is referred to as invariant or rotation invariant if 

P(V- 1 HV) = P(H) (1.1.1) 
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where V is a fixed element in the group diagonalizing H, that is, in SO(iV), 
SU(iV) or USp(2iV) for (3 = 1,2,4, respectively. Equation ([LLT[) implies that 
the probability density function only depends on the eigenvalues, 

P(H) = P(X) = P(x 1 ,...,x N ). (1.1.2) 

Here, we write the diagonalization of the random matrix as H = U~ 1 XU with 
X = diag (xi, . . . , xjy) for j3 = 1, 2 and X = diag (xi,xi, . . . ,xn,xn) for /3 = 
4. The fc-point correlation function R k (xi, . . . ,x k ) measures the probability 
density of finding a level around each of the positions x±, . . . , x k , the remaining 
levels not being observed. One has |Dys62 [Meh04j 



R k ( Xl ,...,x k )= (N m k)l f dx k+1 --- f dx N \A N (X)\Pp(X) , (1.1.3) 

— oo — oo 

where Atv(^) is the Vandermonde determinant. If the probability density 
function factorizes, 

N 

^ = 11^). (1-1-4) 

n=X 

with a probability density function P^ E \x n ) for each of the eigenvalues, the 
correlation functions ([1.1. 3p can be evaluated with the Mehta-Mahoux theo- 
rem |Meh04j . They are k x k determinants for f3 = 2 and 2k x 2k quaternion 
determinants for /3 = 1,4 whose entries, the kernels, depend on only two of the 
eigenvalues x\, ■ ■ ■ ,x k . 

Formula ([1.1. 3p cannot serve as the starting point for the Supersymme- 
try method. A reformulation employing determinants is called for, because 
these can be expressed as Gaussian integrals over commuting or anticommut- 
ing variables, respectively. The key object is the resolvent, that is, the ma- 
trix {x~ — H)' 1 where the argument is given a small imaginary increment, 
Xp = x p — ie. The fe-point correlation functions are then defined as the en- 
semble averaged imaginary parts of the traces of the resolvents at arguments 
x\ , . . . , x k} 

k 

R k (x 1 ,...,x k ) = ^ J : f P(H)f[lmtv _ 1 d[H) . (1.1.5) 

7T K J x *- x p - H 

The necessary limit e — > is suppressed throughout in our notation. We write 
d[-] for the volume element of the quantity in square brackets, that is, for the 
product of the differentials of all independent variables. The definitions ([1.1. 3p 
and (|1.1.5p are equivalent, but not fully identical. Formula ([1.1. 5p yields a sum 
of terms, only one coincides with the definition ([1.1. 3p . all others contain at 
least one 5 function of the form 5(x p — x q ), see Ref. [Guh98]. 
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Better suited for the Supersymmetry method than the correlation func- 
tions ([1.1. 5p are the correlation functions 

k 

Rk(x u ...,x k ) = -\ fd[H]P{H)f[tT -J (1.1.6) 

TV J _ Xp l±jp£ tl 

which also contain the real parts of the resolvents. The correlation func- 
tions (J1.1.5P can always be reconstructed, but the way how this is conveniently 
done differs for different variants of the Supersymmetry method. In Eq. (|1.1.5|) . 
all imaginary increments are on the same side of the real axis. In Eq. (|1.1.6|) . 
however, we introduced quantities L p which determine the side of the real axis 
where the imaginary increment is placed. They are either +1 or —1 and define 
a metric L. Hence, depending on L, there is an overall sign in Eq. ([1.1. 6|) which 
we suppress. We use the short hand notations x^ = x p — iL p e in the sequel. 
In some variants of the Supersymmetry method, it is not important where the 
imaginary increments are, in the supersymmetric non-linear a model, however, 
it is of crucial importance. We return to this point. 

To prepare the application of Supersymmetry, one expresses the correlation 
functions ([1.1. 6p as derivatives 



1 d k 

R k ( Xl , ...,x k )= ( 27r )fc T-rfc qj Zk{ - X + J ) 



(1.1.7) 



of the generating function 

*<* + .,) - jm wn c:g::::":i;g ) 7 ^ 

with respect to source variables J p , p = 1, . . . , k. For j3 = 1, 2 one has 7 = 1 
whereas 7 = 2 for f3 = 4. For later purposes, we introduce the 2k x 2k matrices 
x = diag(xi,xi, . . . ,x k ,Xk) and J = diag(— Ji,+Ji, . . . ,-Jk,+Jk) for /3 = 2 as 
well as the Ak x 4A; matrices x = diag (x±, x\, x\, x\ . . . , x k ,x k ,x k ,x k ) and J = 
di&g(-Ji,-Ji,+Ji,+Ji,. . . ,-J k ,-Jk,+Jk,+Jk) for ^ = 1,4, which appear 
in the argument of Zj~. We write x = x — iLe. Importantly, the generating 
function is normalized at J = 0, that is, Zk(x) = 1. 

1.2 Supermathematics 

Martin [Mar59] seems to have written the first paper on anticommuting vari- 
ables in 1959. Two years later, Berezin introduced integrals over anticommut- 
ing variables when studying second quantization. His posthumously published 
book [Ber87j is still the standard reference on supermathematics. 
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1.2.1 Anticommuting Variables 

We introduce Grassmann or anticommuting variables £ p , p = 1, . . . , k by re- 
quiring the relation 

(pCq = -CqCp , p,q = l,...,k . (1-2.1) 

In particular, this implies ( p = 0. These variables are purely formal objects. In 
contrast to commuting variables, they do not have a representation as numbers. 
The inverse of an anticommuting variable cannot be introduced in a meaningful 
way. Commuting and anticommuting variables commute. The product of an 
even number of anticommuting variables is commuting, 

VSpCqJsr = S>p(,(j(,r = — (,p(,r(,<j = +SrSpS</ = Sr \spsq) • [L.Z.Z) 

We view the anticommuting variables as complex and define a complex conjuga- 
tion, C is the complex conjugate of C, p . The variables Q p and C are independent 
in the same sense in which an ordinary complex variable and its conjugate are 
independent. The property (jl.2.ip also holds for the complex conjugates as well 
as for mixtures, ( P Q = —Q( p - There are two different but equivalent ways to 
interpret (C»)*- The usual choice in physics is 

«;y = c* = -c P , P = i,...,k, (i.2.3) 

which has to be supplemented by the rule 

(C P Cr--Cr)* = c;C g * •••£• (1-2-4) 

There is a concept of reality, since we have 

(c;c P r = c P %* = -CpC P = c;c P • (1.2.5) 

Hence, we may interpret ( p Cp as the modulus squared of the complex anti- 
commuting variable £ p . Alternatively, one can use the plus sign in Eq. (11. 2. 3ft 
and reverse the order of the anticommuting variables on the right hand side of 
Eq. (|1.2.4p . In particular, this also preserves the property (jl.2.5|) . 

Because of C p = and since inverse anticommuting variables do not exist, 
functions of anticommuting variables can only be finite polynomials, 

/(Ci,...,a,cr,...,a)= £ f mi ... mkh ...i k cr---cr(Ci) h ---(ay h (1.2.6) 

m B =0,l 



with commuting coefficients fmi—m h h—l h - Thus, just like functions of matrices, 
functions of anticommuting variables are power series. For example, we have 

exp «C P ) = 1 + aC p C P = ,* - (1-2.7) 

1 "SpSp 

where a is a commuting variable. 
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1.2.2 Vectors and Matrices 

A supermatrix a is defined via block construction, 



a 



a n 
v b 



1.2. 



where a and b are matrices with ordinary complex commuting entries while the 
matrices \x and v have complex anticommuting entries. Apart from the restric- 
tion that the blocks must match, all dimensions of the matrices are possible. 
Of particular interest are quadratic k\jki x k\/k2 supermatrices, that is, a and 
b have dimensions k\ x k\ and k<i x /c2, respectively, /i and v have dimensions 
k\ x /c2 and k% X k\. A quadratic supermatrix a can have an inverse a . 
Equally important are supervectors, which are defined as special supermatrices 
consisting of only one column. As seen in Eq. (jl.2.8|) . there are two possibilities 



t 



and 



$ 



c 

z 



1.2.9) 



where z is a k\ component vector of ordinary complex commuting entries z p , 
and C is a ^2 component vector of complex anticommuting entries C, p . In the 
sequel we work with the first possibility, but everything to be said is valid 
for the second one accordingly. The standard rules of matrix addition and 
multiplication apply, if everything in Sec. I1.2.T1 is taken into account. Consider 
for example the supervector if)' given by 



ijj = (Tip 



a fi 
v b 



z 

c 



az + fi( 
vz + bC, 



(1.2.10) 



which has the same form as ip. Hence the linear map (11.2. 10j) transforms com- 
muting into anticommuting degrees of freedom and vice versa. 

The transpose a and the Hermitean conjugate a> are defined as 



a 



/" 



—v 
b T 



and 



a' 



T\* 



[v 1 ) 



(1.2.11) 



(J2&1 carries over to 



The minus sign in front of v T ensures that (o\(J2) r 
supermatrices o\ and 02- Importantly, {cry = o always holds, but (a T ) T is in 
general not equal to a. As a special application, we define the scalar product 
ijj^X where each of the supervectors ip and x nas either the first or the second 
of the forms (|1.2.9p . Because of the reality property (|1.2.3j) . the scalar product 
ip'ip is real and can be viewed as the length squared of the supervector ifi. 
To have cyclic invariance, the supertrace is defined as 



str a = tr a — tr b 



(1.2.12) 



6 
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such that strait = stromal for two different supermatrices o\ and oi- Corre- 
spondingly, the superdeterminant is multiplicative owing to the definition 



sdet a 



det (a — \ib l i 



det a 



det b det (b — va 1 fi) 

for det b ^ such that sdet g\Oi = sdet uisdet a 2 - 



(1.2.13) 



1.2.3 Groups and Symmetric Spaces 

For and introduction to this topic see chapter 3. Here we only present the 
salient features in the context of the supersymmetry method. The theory of 
Lie superalgebras was pioneered by Kac [Kac77| . Although the notion of super- 
groups, particularly Lie supergroups, seems to be debated in mathematics, a 
consistent definition from a physics viewpoint is possible and — as will become 
clear later on — urgently called for. All supermatrices u which leave the length 
of the supervector ip invariant form the unitary supergroup U(&i/&2). With 
ip' = uip we require {ip'yij/ = ^u^uip = ^ip and the corresponding equation 
for tp' = v)ijj. Hence we conclude 



iru 



1 



nil 



1 



and thus 



u 



u 



-i 



(1.2.14) 



The direct product U(fci) X U(&2) of ordinary unitary groups is a trivial sub- 
group of U(ki/k2), found by simply putting all anticommuting variables in u to 
zero. Non-trivial subgroups of the unitary supergroup exist as well. Consider 
commuting variables, real w p , p = 1, . . . , k\ and complex z p j, p = 1, . . . , k±, j = 
1,2. We introduce the real and quaternion-real supervectors 



W\ 



i, 



Wkx 

Ci 



Cfc 2 



c 



k 2 



and 



$ 



Z\\ 


-*h 


Zl2 


4i 


Zfe x l 




Zk ± 2 


* 

z fcii 


Cl 


-Ci 



a 



L ^fe 2 



"Cfc 2 



(1.2.15) 



The unitary-ortho-symplectic subgroup of the unitary supergroup leaves the 
lengths of ip invariant: UOSp(A;i/2A;2) the length of the first and UOSp(2fci/A;2) 
the length of the second supervector in Eq. (jl.2.15p . Due to the quaternion 
structure in the commuting entries of the second supervector, the proper scalar 
product reads txip^ip. The trivial ordinary subgroups are 0{k\) x USp(2A;2) 
C UOSp(A:i/2A;2) and USp(2A:i) x 0{k 2 ) C UOSp(2A:i/A:2), respectively. 
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As in the ordinary case, non-compact supergroups result from the require- 
ment that the bilinear form ijyLij) remains invariant. The metric L is without 
loss of generality diagonal and only contains ±1. We then have u^Lu = L. 

A Hermitean supermatrix a is diagonalized by a supermatrix u G ~U(ki/k 2 ), 

a = u' 1 su with s = diag(su, . . . ,s kll ,si 2 , . . . , Sfc 22 ) . (1.2.16) 

All eigenvalues s p j are real commuting. Zirnbauer |Zir96aj gave a classification 
of the Riemannian symmetric superspaces. The Hermitean symmetric super- 
space is denoted A| A. Of interest are also the symmetric superspaces AI| All and 
AII|AI. The former consists of the k\j2k 2 x k\/2k 2 supermatrices a = u~ 1 su 
with u G UOSp(A;i/2fc2) and with s = diag (sn, . . . , s kll ,si 2 , s± 2 , • • • , s k22 , s k22 ), 
the latter of the 2k\/k 2 X 2k\/k 2 supermatrices a with u G UOSp(2A;i/A:2) and 
with s = diag (su, sn, . . . , s kl i, s kl i, si 2 , . . . , s k22 ). 

1.2.4 Derivatives and Integrals 

Since anticommuting variables cannot be represented by numbers, there is noth- 
ing like a Riemannian integral over anticommuting variables either. The Berezin 
integral [Bcr87] is formally defined by 

f d( p = and J C P dC P = -^= , (1.2.17) 

and accordingly for the complex conjugates £*. The normalization involving 
v27r is a common, but not the only convention used. The differentials dC, p have 
all the properties of anticommuting variables collected in Sec. 11.2.11 Thus, the 
Berezin integral of the function ()1.2.6[) is essentially the highest order coefficient, 
more precisely /i...n...i/(27r) apart from an overall sign determined by the 
chosen order of integration. For example, we have 

II exp (aC p *C P ) d( p dC p = ^ ■ (1-2.18) 

This innocent-looking formula is at the heart of the Supersymmetry method: 
Anticipating the later discussion, we notice that we would have found the inverse 
of the right hand side for commuting integration variables z p instead of £ p . 

One can also define a derivative as the discrete operation d( p /dC q = 5 pq . To 
avoid ambiguities with signs, one should distinguish left and right derivatives. 
Obviously, derivative and integral coincide apart from factors. Mathematicians 
often prefer to think of the Berezin integral as a derivation. In physics, however, 
the interpretation as integral is highly useful as seen when changing variables. 
We first consider the k 2 vectors £ and r/ = aC, of anticommuting variables 
where a is an ordinary complex k 2 x k 2 matrix. From the definition (J1.2.17P 
we conclude d[rj\ = det _1 ad[£]. This makes it plausible that the change of 
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variables in ordinary space generalizes for an arbitrary transformation \ = x(V') 
of supervectors in the following manner: Let y be the vector of commuting and 
rj be the vector of anticommuting variables in x> then we have 



Qx 
d[x] = sdet — y d[ip\ = sdet 



dy/dz T dy/dC, T 
dr//dz T dr]/d( T 



d[tp] (1.2.19) 



with d[x] = d[y]d[r]] and d[ip] = d[z]d[(}. The Jacobian in superspace is referred 
to as Berezinian. Absolute value signs are not needed if we agree to only trans- 
form right-handed into right-handed coordinate systems. Changes of variables 
in superspace can lead to boundary contributions which have no analog in or- 
dinary analysis. In physics, they are referred to as Efetov-Wegner terms, see 
Rcf. |Rot87j for a mathematical discussion. 

Importantly, the concept of the 5 function has a meaningful generalization 
in superspace. An anticommuting variable £ p acts formally as 5 function when 
integrating it with any function f(( p ), hence 5(C P ) = V2tt( p . More complicated 
are expressions of the form 5(y — <^C) with an ordinary commuting variable 
y and a k component vector of complex anticommuting variables C,- To make 
sense out of it, it has to be interpreted as 

k 

S(y - C f C) = E { —r^ siK) (y)(Ckr • (1-2.20) 

This is a terminating power series, because (C^C) K = for k > k. 

1.3 Supersymmetric Representation 

Several problems in particle physics would be solved if each Boson had a 
Fermionic and each Fermion had a Bosonic partner. A review of this Super- 
symmetry can be found in Rcf. [Mar05 \ Although mathematically the same, 
Supersymmetry in condensed matter physics and Random Matrix Theory has 
a completely different interpretation: the commuting and anticommuting vari- 
ables do not represent Bosons or Fermions, that is, physical particles. Rather, 
they are highly convenient bookkeeping devices making it possible to drastically 
reduce the number of degrees of freedom in the statistical model. Since as many 
commuting as anticommuting variables are involved, one refers to it as Super- 
"symmetry" — purely formally just like in particle physics. In 1979, Parisi and 
Sourlas [Par79j introduced superspace concepts to condensed matter physics. 
Three years later, Efetov [Efe82] constructed the supersymmetric non-linear a 
model for the field theory describing electron transport in disordered systems. 
Efetov and his coworkers developed many of the tools and contributed a large 
body of work on Supersymmetry [Efe83]. The first applications of Supersym- 
metry to random matrices, that is, in the language of condensed matter physics, 
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to the zero- dimensional limit of a field theory, were given by Verbaarschot and 
Zirnbauer [Ver85a] and by Verbaarschot, Zirnbauer and Weidenmuller |Ver85bj . 
Reviews can be found in Refs. |Efe97[ IGuh98[ IMirOOj . see also the chapters on 
chiral Random Matrix Theory and on scattering. 

1.3.1 Ensemble Average 

Using Supersymmetry, the ensemble average in the generating function ([1.1. 8p 
is straightforward. We begin with the unitary case (3 = 2 and express the 
determinants as Gaussian integrals 

— j — = / d[z p ] exp (iLpZ^H - x p + J p )z p 



det(H — x p + J p 

'p 



det(H — x^ + J p 



/ d[z p ] exp (iLpZ^H 
J d[C p ] exp (i$(H -x±- J p )Cp) (1.3.1) 



(2tt) n 

over altogether k vectors z p , p = 1, . . . , k with TV" complex commuting entries 
and k vectors £ p , p = 1, . . . , k with N complex anticommuting entries. When 
integrating over the commuting variables, the imaginary increment is needed 
for convergence, for the integrals over anticommuting variables, convergence 
is never a problem. Hence we may write the metric tensor in the form L = 
diag (Li, . . . , Lfc, 1, . . . , 1). Collecting all H dependences, the ensemble average 
in Eq. (|1.1.8j) amounts to calculating 

*(#)= J d[H]P(H)exp(itiHK) . (1.3.2) 

where the N x N matrix K assembles dyadic products of the vectors z p and C 



v 



k 



K = y(LpZp4-(pCl) . (1.3.3) 



/ j y-'p'Wp spSp 
P =i 

For all L, this is a Hermitean matrix K^ = K. 

We now turn to the orthogonal case f3 = 1. At first sight it seems irrelevant 
whether H is Hermitean or real-symmetric in the previous steps. However, the 
Fourier transform ([1.3.21) only affects the real part of K, because the imaginary 
part of K drops out in trHK if H is real-symmetric. Thus, instead of the 
Gaussian integrals (|1.3. lj) , we rather use 

r N r 

d[w ( V] exp (iL pWp ^ T {H -x± + J p ) W V> 



■n 



det(H -xp + J p ) 



/' 



d[w p '] exp ( iLpiu^ ' (H — x p + J p )w. 



det{H - xf + Jp) f (-Aitj ± 



n N 



d[C P ] exp [iCl{H -Xp- J p )( p j 

e W (-i^(H-x^-J p )C p ) , (1-3.4) 



10 
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(i) 



,(2) 



where the to N component vectors w p and w p have real entries. For each 
p, we can construct a 4N component supervector out of w p , w p , ( p and C 
whose structure resembles the one of the first of the supervectors (|1.2.15|) . but 
with a different number of components. Reordering terms, we arrive at the 
Fourier transform f 1 1 . 3 . 2 [) . but now for real-symmetric H and with 



K = ]T [L p w^w^ T + L p w^w^ T - C^ + C P C P ) , 



1.3.5) 



which is N x N real-symmetric as well. For f3 = 4, one has to reformulate the 
steps in such a way that the corresponding K becomes self-dual Hermitean. 

1.3.2 Hubbard-Stratonovich Transformation 

Due to universality, it suffices to assume a Gaussian probability density function 
P{H) ~ exp(— /3tr H 2 /2) in almost all applications in condensed matter and 
many-body physics as well as in quantum chaos. Hence the random matrices 
are drawn from the Gaussian orthogonal (GOE), unitary (GUE) or symplectic 
ensemble (GSE). The Fourier transform (J1.3.2J) is then elementary and yields a 
Gaussian. The crucial property 



<&(K) = exp 



■k lK ' 



exp 



-k ixB '' 



1.3.6) 



holds, where B is supermatrix containing all scalar products of the vectors to 
be integrated over. The second equality sign has a purely algebraic origin. For 
f3 = 2, B has dimension k/k x k/k and reads 



B = L 1 ' 2 



\z\ 


z\ Z k 


4(i ■ ■ 


z\Ck 


i z i ■ ■ 


■ z l z k 


4ci •• 


■ 4<k 


chi ■■ 


■ -C\zk 


-clci •• 


■ -ClCk 


cUi ■■ 


■ -c,Ufc 


-clci •• 


■ -cjfc* 



L l/2 



1.3.7) 



While K is Hermitean, the square roots L 1 ' 2 destroy this property for B, since 
L 1 2 can have imaginary units % as entries, B is Hermitean only for L = 1. In 
general, B is in a deformed (non-compact) form of the symmetric superspace 
A|A. For /3 = 1,4, the supermatrix B has dimension 2k/2k x 2k /2k and it 
is in deformed (non-compact) forms of the symmetric superspaces AI|AII and 
AII|AI, respectively. We give the explicit forms later on. The identity f 1 1 . 3 . 6 [) 
states the keystone of the Supersymmetry method. The original model in the 
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space of ordinary N x N matrices is mapped onto a model in space of superma- 
trices whose dimension is proportional to k, which is the number of arguments 
in the fc-point correlation function. 

The Gaussians (11.3.6P contain the vectors, that is, their building blocks to 
fourth order. To make analytical progress, a Hubbard-Stratonovich transfor- 
mation in superspace is used, 

exp f-i-strS 2 ^ = C W /exp (-|str (La)A exp (istr l}l 2 cL x l 2 B\ d[a] , 

(1.3.8) 
where c^ = 2 k ( k ~^ and c^ = 2 k ( Ak ~ 3 ^ 2 for @ = 1,4. We notice the appear- 
ance of the matrices L and L 1 ' 2 in (jl.3.8p . For L = 1, the supermatrices a 
and B have the same symmetries. However, as already observed in the early 
eighties for models in ordinary space [Sch80, Pru82], this choice is impossible 
for L/l, because it would render the integrals divergent. There are two ways 
out of this problem. One either constructs a proper explicit parameterization 
of a or one inserts the matrices L and L 1 ' 2 according to (jl.3.8p . A mathemat- 
ically satisfactory understanding of these issues was put forward only recently 
in Ref. |Fyo08| . 

Another important remark is called for. Because of the minus sign in the su- 
pertrace ()1.2.12|) . a Wick rotation is needed to make the integral convergent. It 
formally amounts to replacing the lower right block of a, that is, b in Eq. (|l,2,8p . 
with ib. Apart from that, the metric L is also needed for convergence reason. 
Now the vectors appear in second order. They can be ordered in one large 
supervector ^f. For (3 = 2 it has the form (J1.2.9P with k\ = k^ = kN, for /? = 1 
it has the first of the forms (J1.2.15H with ki = 2kN, k 2 = kN and for /3 = 4 it 
has the second of the forms (|1.2.15|) with k\ = kN, k 2 = 2kN. The integral to 
be done is then seen to be the Gaussian integral in superspace 

/exp (Wt (l}l^{l}l 2 aL x l 2 -x ± - J)L 1 ' 2 ® 1^) tf) d[tf] 

= sdet ~ N ^ 2 ^(aL -x ± -J), (1.3.9) 
where the power A^ is due to the direct product structure. We eventually find 

Z k (x + J) = c^ f ^exp T-^str {La) 2 ) sdet~ N ^ 2 ^ (aL - x ± - J)d[a] (1.3.10) 

as super symmetric representation of the generating function. The average over 
the N x N ordinary matrix H has been traded for an average over the matrix 
a whose dimension is proportional to k, that is, independent of N. 

1.3.3 Matrix 5 Functions and an Alternative Representation 

In Refs. [Leh951 l"Hac95| . a route alternative to the one outlined in Sec. 11.3.21 was 
taken. These authors used matrix 5 functions in superspace and their Fourier 
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representation to express functions f{B) of the supermatrix B in the form 
f(B)= [f(p)5(p-B)d[p] = c^ 2 fd[p}f(p) fd[<r]exp(-istva(p-B)) , 



(1.3.11) 
where auxiliary integrals over supermatrices p and a are introduced. For sim- 
plicity, we only consider L = 1 here. The function 5{p — B) is the product of 
the 5 functions of all independent variables. As discussed in Sec. 11.2.41 it is 
well-defined. For all functions /, formula (|1.3.1ip renders the integration over 
the supervector $ Gaussian. When studying Gaussian averages of ratios of 
characteristic polynomials, Fyodorov |Fyo02| built upon such insights to con- 
struct an alternative representation for the generating function. He employs a 
standard Hubbard-Stratonovich transformation for the lower right block of the 
supermatrix B in Eq. (jl.3.7|) which contains the scalar products Q>Cq- He then 
inserts a 5 function in the space of ordinary matrices to carry out the integrals 
over the vectors z p . Although Supersymmetry is used, the generating function 
is finally written as an integral over two ordinary matrices with commuting 
entries. In this derivation, the Ingham-Siegel integral 

/(«*>(«) = J e W (-tTRS)det m Sd[S] ~ ^j, R (1.3.12) 

s>o 
for ordinary Hermitean N x N matrices R and S appears, where m > 0. 

1.3.4 Generalized Hubbard-Stratonovich Transformation and 
Superbosonization 

Is Supersymmetry only applicable to Gaussian probability density functions 
P(H) ? - - In Ref. [Hac95j, Supersymmetry and asymptotic expansions were 
used to prove universality for arbitrary P{H). The concept of superbosoniza- 
tion was put forward in Ref. [Efe04] and applied in Ref. [Bun07| to a gener- 
alized Gaussian model comprising a variety of correlations between the ma- 
trix elements. Extending the concept of superbosonization, a full answer to 
the question posed above was given in two different but related approaches in 
Refs. |Guh06, Kie09a] and [Lit08]: An exact super symmetric representation ex- 
ists for arbitrary, well-behaved P(H). As the equivalence of the two approaches 
was proven in Ref. [Kie09b], we follow the line of arguing in Refs. [Guh06, 
IKie09a| . For (3 = 2, we define the N x 2k rectangular supermatrix 

A=[zf-z k Ci---Ck] , (1-3.13) 

where the z p , p = 1, . . . , k and £„, p = 1, . . . , k are A^ component vectors with 
complex commuting and anticommuting entries, respectively. We also define 
the N x 4k supermatrix 

A=[z 1 zl---z k zlCiCi--<kC k ] (1-3.14) 
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for p = 1 and eventually the 2N x Ak supermatrix 



A 



-4 


Zk 


* 

-4 


Ci 


-Ci* 


a 


-a 


4 


Zk 


z* 

z k 


Ci 


cr 


a 


c k 



(1.3.15) 



for (3 = 4. This enables us to write the ordinary matrix K introduced in 
Sec. II. 3. H and the supermatrix B introduced in Sec. 11.3.21 for all f3 in the form 

K = ALA^ = {AL 1/2 ){L 1/2 A ] ) 

B = (L l l 2 A^){AL l l 2 ) = L l l 2 A^AL l l 2 . (1.3.16) 

For /3 = 2, we recover Eq. (11.3.7p . This algebraic duality between ordinary and 
superspace has far-reaching consequences. One realizes [Guh91[ IGuh06[ ILit08] 
that the integral (II. 3. 2D is the Fourier transform in matrix space of every, ar- 
bitrary probability density function P{H) and that &(K) is the correspond- 
ing characteristic function. Since we assume that P(H) is rotation invariant, 
the same must hold for &(K). Hence, &(K) only depends on the invariants 
tr K m , m = 1, 2, 3, ... . Due to cyclic invariance of the trace, the duality (jl .3.16|) 
implies for all m the crucial identity 

trK m = strB m , such that $(K) = $(B) . (1.3.17) 

Hence, viewed as a function of the matrix invariants, $ is a function in ordinary 
and in superspace. We now employ formula ([1.3. lip for &(K) = &(B), do the 
Gaussian ^ integrals as usual find for the generating function 

Z k (x + J) = c (/3)2 I exp (-itr (x + J)Lp) ${p)I(p)d[p] (1.3.18) 

with I(p) being a supersymmetric version of the Ingham-Siegel integral. The 
supermatrices p and a have the same sizes and symmetries as B. A convolution 
theorem in superspace yields the second form 

Z k (x + J)= J n(a)sdet ~^/ 2 7 (aL - x ± - J) d[a] , (1.3.19) 

where n(<j) is the superspace Fourier backtransform of the characteristic func- 
tion <£(/>). It plays the role of the probability density function for the super- 
symmetric representation. To apply these general results for exact calculations, 
explicit knowledge of either $(/?) or n(<r) is necessary. 

1.3.5 More Complicated Models 

Most advantageously, Supersymmetry allows one to make progress in impor- 
tant and technically challenging problems beyond the invariant and factorizing 
ensembles, for example: 
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• Invariant, but non-f actor izing ensembles. The probability density func- 
tion P(H) has the property ([1.1. ip . but not the property (|1.1.4|) . They 
can, in principle, be treated with the results of Sec. 11.3.41 



Sparse or banded random matrices |Fyo91 [Mir91| . see chapter 23. The 



probability density function P(H) lacks the invariance property (jl.l.ip . 

• Crossover transitions or external field models. One is interested in the 
eigenvalue correlations of the matrix H(a) = H^°> + aH, where H is a 
random matrix as before and where H^' is either a random matrix with 
symmetries different from H or a fixed matrix. The parameter a measures 
the relative strength. As the resolvent in question is now (x~ — H(a)) = 
(Xp — (H(°> +aH))~ l , we have to replace H by H(a) in the determinants 
in Eq. (|1.1.8|) . but not in the probability density function P(H) which 
usually is chosen invariant, see Ref. [Guh96bJ. 

• Scattering theory and other problems, where matrix elements of the resol- 
vents enter |Ver85b| , see chapter 2 on history and chapter 34 on scattering. 
In the Heidelberg formalism [Mah69], scattering is modeled by coupling 
an effective Hamiltonian which describes the interaction zone to the scat- 
tering channels. The resolvent is then (x~ — H + iW)~ l where the N x N 
matrix W contains information about the channels. One has to calculate 
averages of products of matrix elements \{x~ — H + iW / )~ 1 ] nm . To make 
that feasible, the source variables J p have to be replaced by N x N source 
matrices J p and instead of the derivatives (jl.l.7|) . one must calculate 
derivatives with respect to the matrix elements J p>nm . The probability 
density function P(H) is unchanged. 

• Field theories for disordered systems, see Refs. |Efe83[lEfe97| . 

Of course, these and other non-invariant problems can not only be studied 
with the Supersymmetry method, other techniques ranging from perturbative 
expansions, asymptotic analysis to orthogonal polynomials supplemented with 
group integrals are applied as well, see chapters 4., 5. and 6. Nevertheless, the 
drastic reduction in the numbers of degrees of freedom, which is the key feature 
of Supersymmetry, often yields precious structural insights into the problem. 

1.4 Evaluation and Structural Insights 

To evaluate the supersymmetric representation, a large iV expansion, the cel- 
ebrated non-linear a model, is used in the vast majority of applications. We 
also sketch a method of exact evaluation which amounts to a diffusion process 
in superspace. Throughout, we focus on the structural aspects. A survey of the 
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numerous results for specific systems is beyond the scope of this contribution, 
we refer the reader to the reviews in Refs. |Efc97, Guh98, MirOO]. 

1.4.1 Non-linear a Model 

The reduction in the numbers of degrees of freedom is borne out in the fact that 
the dimension N of the original random matrix H is an explicit parameter in 
Eqs. (|1.3.1U|) and (|1.3.19|) . Hence we can obtain an asymptotic expansion in 1/N 
by means of a saddle point approximation [Efe83[ Efe9~7" l IVer85al IVer85b| . This 
suffices because one usually is interested in the correlations on the local scale 
of the mean level spacing. Hence, the saddle point approximation goes hand in 
hand with the unfolding. The result of this procedure is the supersymmetric 
non-linear a model. We consider the two-point function k = 2. The integrand 
in Eq. (J1.3.10J) is written as exp(— F(x + </)) with the free energy 

F(x + J) = ^str (La) 2 + -/str \n{oL - x ± - J) , (1.4.1) 

2 27 

which is also referred to as Lagrangean. We introduce center x = (x\ + #2)/2 
and difference Ax = X2 — x\ of the arguments. In the large N limit, £ = Ax/D 
has to be held fixed where D ~ 1/y/N is the local mean level spacing. Hence, 
when determining the saddle points, we may set Ax = such that x = xl. 
Moreover, as we may choose the source variables arbitrarily small, we set J = 
as well. Since all symmetry breaking terms are gone, the free energy F(x) 
with x = xl is invariant under rotations of a which obey the metric L. Thus, 
variation of F[x) with respect to a yields the scalar equation 



sq{x — sq) = — , such that so = — I x ± i\\ x 2 ) (1-4-2) 

27 2 \ 




inside the spectrum, \x\ < y/2N/j. This is the famous Pastur equation and its 
solution so |Pas72j . The latter is proportional to the large iV one-point func- 
tion whose imaginary part is the Wigner semicircle. We arrive at the important 
insight that the one-point function provides the stable points of the supersym- 
metric representation, the correlations on the local scale are the fluctuations 
around it. To make this more precise, we recall the result of Schafer and Weg- 
ner |Sch80| for the non-linear a model in ordinary space: when doing the large 
N limit as sketched above, the imaginary increments of the arguments x\ and 
X2 must lie on different sides of the real axis. Otherwise, the connected part 
of the two-point function cannot be obtained as seen from a contour-integral 
argument. Hence, the metric L must not be proportional to the unit matrix, 
the groups involved are non-compact and a hyperbolic symmetry is present. 
This carries over to the commuting degrees of freedom in superspace [Efe82j . 
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the groups are U(l,l/2) for /3 = 2 and UOSp(2,2/4) for /3 = 1,4. The full 
saddle point manifold is found to be given by all non-compact rotations of 
do = x/2 + iy / 2N/^f — x 2 L/2 which leave F(xl) invariant. One parameterizes 
the group as u = uov with uq in the direct product U(l/1) x U(l/1) for /3 = 1 
and in UOSp(2/2) x UOSp(2/2) for /3 = 1, 4 and with v in the coset 



U(l,l/2) 



U(l/1) x U(l/1) 
UOSp(2,2/4) 
UOSp(2/2) x UOSp(2/2) 



for /3 = 2 , 

for ,0 = 1,4. (1.4.3) 



As uq and L commute, the saddle point manifold is v ~ 1 ctqv , that is, essentially 
Q = v Lv with the crucial property Q 2 = 1. To calculate the correlations, 
we must re-insert the symmetry breaking terms Ax and J into the free energy. 
We put a = v~ 1 aov + 5a and expand to second order in the variables 5a which 
are referred to as massive modes. They are integrated out in the generating 
function (jl.3.10p as Gaussian integrals. One is left with integrals over the coset 
manifold, that is, over the Goldstone modes. On the unfolded scale, the two- 
point correlation functions (J1.1.3J) acquire the form 1 — 12(C)- The two-level 
cluster functions read 

y 2 (£) = -Re / exp(^strQL) strM 1 QLstrM 2 QLdn(Q) , (1.4.4) 

where d[i{Q) is the invariant measure on the saddle point manifold. The matri- 
ces Mi, i = 1, 2 result from the derivatives with respect to the source variables, 
Mi is found by formally setting Jj = 1 and Ji = 0, I 7^ i, in the matrix J. The 
expressions (jl.4.4p can be reduced to two radial integrals on the coset mani- 
fold for the GUE and to three such integrals for GOE and GSE. Efetov |Efe83| 
discovered Eq. (jl.4.4p when taking the zero-dimensional limit of his super- 
symmetric non-linear a model for electron transport in disordered mesoscopic 
systems. He thereby established a most fruitful link between Random Matrix 
Theory and mesoscopic physics. 

The non-linear a model, particularly its mathematical aspects, was recently 
reviewed in Ref. |Zir06| . 

1.4.2 Eigenvalues and Diffusion in Superspace 

The supersymmetric representation and Fyodorov's alternative representation 
of Sec. 11.31 are exact for finite N. In some situations, it is indeed possible 
and advantageous to evaluate them without using the non-linear a model. It 
has been shown for the supersymmetric representation [Guh 06, Kic09a that 
the imaginary increments of the arguments may then all lie on the same side 
of the real axis. We have L = 1 and all groups are compact. As we aim 
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at structural aspects, we consider the crossover transitions involving H{a) = 
H^> +aH as discussed in Sec. 11.3,51 We introduce the fictitious time t = a> 2 /2. 
Dyson |Dys62, Dys72| showed that the eigenvalues of H{t) follow a Brown- 



ian motion in t which implies that their probability density function is prop- 
agated by a diffusion equation. Without any loss of information, Supersym- 
metry reduces this stochastic process to a Brownian motion in a much smaller 
space [Guh96bj. It is precisely the radial part of the Riemannian symmet- 
ric superspaces discussed in Sec. 11.2.31 The quantity propagated is then the 
generating function Z k {x + J,t) of the correlations. The initial condition 

Zf\s)= f P^ (H^)sdet- 1 (s® 1 + 1® H^)d[H( 1] (1.4.5) 

is arbitrary, it includes ensembles, but also a fixed matrix H^°> if the probability 
density P(°\H(°>) is chosen accordingly. Due to the direct product structure, 
Zi (s) is rotation invariant. The diagonal matrix s is in the above mentioned 
radial space, such that a = u~ 1 su, see Sec. 11.2.31 For (5 = 2, this space coincides 
with x + J, for (5 = 1,4, it is slightly larger. The generating function is then a 
convolution in the radial space, 

Z k (r,t) = j T k {s,r,t)Z^\s)B k { S )d[s] . (1.4.6) 

When going to eigenvalue-angle coordinates a = u~ l su Berezinians B k (s) occur 
analogous to |Aat(X)| /3 in ordinary space. The propagator is the supergroup 
integral 

T k (s,r,t) = c^' exp I — — str (s + r ) J / exp I — — strn~ sur J d/j,(u) . 

(1.4.7) 
For (5 = 2 and all k, this integral is known explicitly [Guh91, Guh96aJ. Unfor- 
tunately, for (5 = 1, 4 the available result [Guh02j is handy only for k = 1 but 
cumbersome for k = 2. 

It is a remarkable inherent feature of Supersymmetry that the propagator 
and thus the diffusion process of Z k {r,t) on the original scale in t carries over 
unchanged to the unfolded scale when introducing the proper time r = t/D 2 . 
The initial condition is the unfolded large N limit of Z\ (s). Moreover and in 
contrast to the hierarchical equations for the correlation functions [Fre88], the 
Brownian motion in superspace for the generating functions is diagonal in k. 

1.5 Circular Ensembles and Color— Flavor Transfor- 
mation 

In many physics applications, the random matrices H model Hamiltonians, 
implying that they are either real symmetric, Hermitean or quaternion real. 
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Mathematically speaking, they are in the non-compact forms of the corre- 
sponding symmetric spaces. However, if one aims at modeling scattering, it 
is often useful to work with random unitary matrices S, taken from the com- 
pact forms of these symmetric spaces [Dys62] . These are U(N)/0(N), U(iV) 
and XJ(2N)/Sp(2N), respectively, leading to the circular orthogonal, unitary 
and symplectic ensemble COE, CUE and CSE which are labeled f3 = 1,2,4. 
The phase angles of S play the same role as the eigenvalues of H. Due to the 
compactness, no Gaussian or other confining function is needed and the proba- 
bility density function is just the invariant measure on the symmetric space in 
question. On the local scale of the mean level spacing, the correlations coincide 
with those of the Gaussian ensembles (Dys62, M eh04] . 

Zirnbauer jZir96bj showed how to apply Supersymmetry to the circular 
ensembles. His approach works for all three symmetry classes, but for simplicity 
we only discuss the CUE which consists of the unitary matrices S = U G XJ(N). 
Consider the generating function 

fe+ , , , . \ T T\ k- 



v ft \ [a <tj\ TT det ( X - ewW+p)U) TT det I 1 ~ ^PW-i)^ 
Z k+k _ (tf , i p) = J dm JUL det(1 _ exp(n?+p)[7) 11 det(l-exp(W_,)tft) ' 

(1-5.1) 

where d/j,(U) is the invariant measure on U(iV). To derive the correlation func- 
tions Rk($i, • • • , ~&k) one sets k + = fe_ = k, takes derivatives with respect to the 
variables <p± p and puts certain combinations of variables ip± p and 'Q± p equal. 
The variables t?± p in Eq. (jl.5.ip have small imaginary increments to prevent 
Zk+k-($i 4 1 ) from becoming singular. 

Since the Hubbard-Stratonovich transformation of Sec. 11.3.21 cannot be 
employed to construct a supersymmetric representation of Zk + k_($,<fi), Zirn- 
bauer [Zir96bj developed the color-flavor transformation based on the identity 



I dfi(U) exp (^l* pj U nm *™ pj + Vl* qj U nm *& 

d[A]d[A]sdet N (l - AA) exp ^ n + * pj A pjql ^ ql + *^A g ^*^).5.2) 



Vq'J) 



which transforms an integral over the ordinary group U(AQ into an integral over 
k + /k + x k-/k- rectangular supermatrices A and A parameterizing the coset 
space XJ(k+ + k-/k+ + k-)/XJ(k+/k+) xU(L/L). The integrals depend on the 
supertensor ^/ with components ^\ p j- The indices n,m = 1,...,N label the 
elements of U in ordinary space. The indices pj and ql are superspace indices 
with p = 1, . . . , k+, q = 1, . . . , k- and with j,l = 1,2 labeling the four blocks 
of the supermatrices, see Eq. (jl.2.8p . Summation convention applies. The 
superfields \& are used to express the determinants in Eq. (|1.5.ip as Gaussian 
integrals. After the color-flavor transformation, they are integrated out again. 
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As the name indicates, the color-flavor transformation is naturally suited 
for applications in lattice gauge theories where U is in the color gauge group, 
see Ref. |Nag01| . In Ref. |Wei0 5] . the color-flavor transformation was derived 
for the gauge group SU(iV) relevant in lattice quantum chromodynamics. 

1.6 Concluding Remarks 

Apart from the wealth of results for specific physics systems which had to be 
left out, some important conceptual issues could not be discussed here either 
due to lack of space: we only mentioned applications of the other Rieman- 
nian symmetric superspaccs [Zir96aJ to Andreev scattering and chiral Random 
Matrix Theory. As random matrix approaches are now ubiquitous in physics 
and beyond, one may also expect that the Supersymmetry method spreads out 
accordingly. From a mathematical viewpoint, various aspects deserve further 
clarifying studies, in the present context most noticeably the theory of super- 
groups and harmonic analysis on superspaces. 
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